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Abstract 

We establish the existence of real pole-free solutions to all even members of the 
(— ! ' Painleve I hierarchy. We also obtain asymptotics for those solutions and describe their 

pi |. relevance in the description of critical asymptotic behavior of solutions to the KdV 

equation in the small dispersion limit. This was understood in the case of a generic 
critical point, and we generalize it here to the case of non-generic critical points. 
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1 Introduction and statement of results 



^ ' We first introduce the Painleve I hierarchy, which has the Painleve I equation 

^ qss = s + 6g2 (1.1) 

O 



as its first member. The m-th equation of the hierarchy is of order 2m and is defined 
recursively. It has the form 



m—l 



s + Cm{q)+^tjCj.i{q)=Q, tl,...,t„_i GM, (1.2) 

• 1—1 . 1=1 

X. 

H where Cm is the Lenard-Magri recursion operator defined by 

Co{q) = -4q, (1.3) 

d [Id? d \ 

^-Cfc+ilQ) = ( ^^ -2g— -gj £yt(g), for A; = 0, ...,m- 1. (1.4) 

The constants of integration in (jl.4p are fixed by the requirements £i(0) = • • • = £m(0) = 
0. One could also add a term tmC-m-i in ()1.2p . but this term cancels after a simple 
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transformation of the other variables q,s,ti, . . . , tm-i- The first equations in the hierarchy 
are given by (up to multiphcation by a non-zero constant) 



m = 
m = 1 
m = 2 
m = 3 



m 



4: 



s - 4g = 0, 

qss = s + Qq^, 

qssss = 4s - 40g3 + lOg^ ^ 20qq,, - IQtiq, 

g(6) = 16s + 28qqssss + ^GQsQsss + 42gL - 280{q^qss + qql - g^) 

+ 16t2(6g^ - qss) - 64tig, 
g(8) = 64s + 36gg(^) + lOSg^g^^^ + 228qssqssss - ^^M^qssss 

+ 138gLs - 1512(7<7L " 1848g2g,, - 2016gg,g,,, 
- 2016g^ + ^'iQQq^qss + f>M<dq^ql 
+ 16t3(-40g3 + 10g2 + 20(7gss - q^^^) 

+ Ut2{<oq^ - qss) - 2bQtiq. 

In the above equations we have written q^^' for the j-th derivative of q with respect to 
s. We win can equation ([OD the P{^ equation, and we refer to pi| [38t liOl H3| [26] 
for more information about the first Painleve hierarchy. Given ti,...,tm-i) solutions 
to these equations are meromorphic functions in the complex s-plane, with in general 
an infinite number of poles \4Q\. Heuristic arguments supporting the existence of real 
pole-free solutions to the even members of the hierarchy were already given in [3] in the 
case where ti = . . . = im-i = 0. Their asymptotic behavior was discussed in [37|. A 
particular solution to the second member of the hierarchy appeared in |41l |42] in relation 
with the Gurevich-Pitaevskii special solutions to the KdV equation [271 [39]. This solution 
corresponds to the solution studied in [3] for ti = 0, and has applications in the study of 
ideal incompressible liquids [32l [33] and quantum gravity [15] , see also [25] . 
After the rescalings 

U = -60^/^ • g, X = m-^/^ -s, T = -4 • 60"^/^ • h, 

the second member of the hierarchy P^ becomes 



X = TU 



1 



1 



U-" + -rm + 2UUxx) + wjK^xxxx 



24 



240 



(1.5) 



This equation was studied by Dubrovin in [16], where he conjectured the existence and 
uniqueness of a real solution without poles for real values of X and T. The existence of such 
a solution was proved in [6j together with the asymptotic behavior U{X,T) ~ =p(6|X|)"^'^ 
as X ^ iboo. In addition it is known that U{X, T) is also a solution to the KdV equation 
Ut + UUx + ^Uxxx = [271 [391 El US]. 

As a first result in this paper, we will prove the existence of real pole-free solutions 
for all even members of the hierarchy, we will obtain asymptotics for them, and show that 
they follow, as functions of the tiine variables ti, . . . ,tm-i, the time flows of the KdV 
hierarchy. 



Theorem 1.1 Let m be an even positive integer. There exists a solution q 
to equation /il.2\) which has the properties 



j{s,ti, . . . ,tm-l) 



(i) q is real and has no poles for real values of s,ti, . . . , tm-i,' for s,ti, . . . , im-i ^^ o, 
sufficiently small neighborhood of the real line, q depends analytically on each of its 
variables, 



lik + oL I 1 —.^k = 0, fork = l,...,m-l, (1.6) 



(ii) q satisfies the PDE 

1 d 

2kTTds' 

which is (up to re-scaling) the k-th equation in the KdV hierarchy, 
(Hi) for ti, . . . , tm-i = 0, q has the asymptotic behavior 

1 __ m 

g(s,0,... ,0) = c|s|™+i +0(151 '"+1), as s ^ ±oo, (1.7) 

with 

sgn(s) /2™"i(m + 1)!\ ^^^+1 



(2m + 1)!! 



{U 



Remark 1.2 The above results are not true for m = 1, since it is known that there do 
not exist real pole-free solutions to the Painleve I equation [1]. Also for m > 1 odd, we do 
not expect that our results can be generalized. Parts (i) and (iii) are the essential parts of 
the theorem, part (ii) will follow from rather standard arguments that express the relation 
between the Painleve I hierarchy and the KdV hierarchy. 

Remark 1.3 For ti, . . . ,tm-i S K fixed, our asymptotic analysis can be generalized to 
obtain the asymptotics 

1 1 

q{s,ti,...,trn~i) = c\s\^+^ +0{\s\ ^+1), as s -> iboo, (1.9) 

see Remark 12.41 below. Note that the error term is weaker in the case of non-zero tj's than 
in (fTTl) . 

Remark 1.4 The leading order of the asymptotic behavior of q as s — )• ±00 for fixed 
ti, . . . , tm-i is relatively simple, and can be formally obtained when neglecting all deriva- 
tives of q in the equations of the hierarchy. It was proved in [8] that solutions to the 
Painleve I hierarchy with this asymptotic behavior exist as x — )■ 00 in a sector containing 
the positive real line. If one takes double scaling limits where ti, . . . , tm-i tend to infinity 
simultaneously with s, the situation becomes more complicated. Then the type of asymp- 
totics will depend on the precise scaling of all variables. For example it can be expected 
that the asymptotics for q can be expressed in terms of elliptic ^-functions in some regions 
and, for ttt, > 2, in terms of hyperelliptic ^-functions in other regions. For critical scalings 
of the variables, one can even expect asymptotics for q = qm va. terms of the pole-free 
solutions q2i<lAi- ■ ■ ■, Qm-2 of the lower order equations in the Painleve I hierarchy, and in 
terms of certain solutions to the Painleve II hierarchy, see the discussion in [32] and [3] . 



1.1 Critical behavior for KdV solutions in the small dispersion limit 

Let us consider the KdV equation 

Mi + 6iiii2,. + e^Ui,.^,.^,. = 0, e > 0. (1-10) 

For small e, this is an example of a Hamiltonian perturbation of the Hopf equation ut + 
Quux = 0. Solutions to the Hopf equation exist only for small times, and develop a point 
of gradient catastrophe after a certain time. Indeed, consider for example initial data uq 
which are negative, smooth, tending to rapidly at itcxD and with a single local minimum. 
Then the method of characteristics describes the solution in terms of the initial data: 

and the slope becomes infinite at the critical time 

1 

t, 



max5giR[-6ii^(^)]' 

The point Xc and time tc where Ux blows up, and the value Uc = u{xc-, tc) are also deter- 
mined by the equations 

F{u;x,t):=-x + Qut + fL{u) = Q, (1.11) 

F'(M;x,t)=6t + /l(M) = 0, (1.12) 

F"{u-x,t) = f'l{u) = Q, (1.13) 

where Jl is the inverse of the decreasing part of the initial data uq. For generic initial 
data we have f'l'{uc) 7^ 0, which means that the Hopf solution behaves locally as 

u{x, tc) = Uc — c{x — Xc)^" + 0{x — Xc), as X — )• Xc- (1-14) 

For non-generic initial data however, it can happen that f'l^{uc) = and that 

1 
u{x, tc) = Uc — c{x — Xc) ""+1 + 0{x — Xc), as X — )• Xc, (1-15) 

for any even value of m. This is the case if the initial data are such that 

/f K) = ft\uc) = ...= f'r\^c) = 0, ft'-'H^c) / 0. (1.16) 

It was conjectured by Dubrovin [16l[T71[l8l[19l[20] that the behavior of generic solutions 
to any Hamiltonian perturbation of a hyperbolic equation near the critical point of the 
unperturbed equation is described universally in terms of the pole-free solution to the 
second member (of order 4) of the Painleve I hierarchy: there should be an expansion of 
the form 

u{x, t, e) = Uc+aie^/'^U (a2e-^'\x -Xc- a3{t - t^)), a4e-^/^(t - tc)) +0 (e^'A , (1.17) 



in a double scaling limit where e — )• and simultaneously x — )• Xc, t — )• tc in such a way 
that e~^''^(x — Xc — a^it — tc)) and e~^''^{t — tc) tend to real constants, and where U is 
the pole- free solution to equation (II. 5p . The values of ai, ... ,04 depend on the equation 
and the initial data, but not on x, t, e. In the case of the KdV equation, this was proved 



afterwards [5j for analytic negative initial data with sufficient decay at ±00 and with a 
single local minimum, under the condition that the gradient catastrophe for the Hopf 
equation is generic, i.e. f'/^'{u;Xc,tc) 7^ 0. If this condition is not satisfied, we will prove 
here that the KdV solution is no longer described in terms of the pole-free solution to 
the Pj equation, but in terms of a solution to a higher order equation in the Painleve 
I hierarchy satisfying the properties given in Theorem II. li The order of the equation is 
determined by the number of vanishing derivatives of fi: if we have (|1.16p for m G 2N, 
the pole-free solution to the Pj" equation of order 2m will appear. 

We consider initial data uq{x) in the class of negative functions with only one local 
minimum, and such that uq can be extended to an analytic function no(-z) in a region of 
the form 

S = {z GC:\Imz\ < tan6l|Rez|} U {z G C : |Imz| < a} 

for some < 6 < 7r/2 and cr > 0. In addition we need sufficient decay at infinity, 

uo{x) = O ( I , s > 0, xeS, X -^ 00. (1-18) 

The local minimum is localized at a point xm and we assume that Uq{xm) 7^ and 
uo{xm) = -1. 

Theorem 1.5 Let uq{x) be initial data for the Cauchy problem of the KdV equation sat- 
isfying the conditions described above, and assume that we have il.lb]) for m € 2N. Write 
Uc = u{xc,tc) for the Hopf solution at the point x^ and time tc of gradient catastrophe of 
the Hopf equation. We take a double scaling limit where we let e ^- and at the same 
time X — )■ Xc and t ^^ tc in such a way that, for some to,ti G M, 

X- Xc- Gucit -tc) .. 3{t - tc) ., , „. 

lim , 2m+2 =^0, -lim , 2 — = n, (1.19) 

^l/2g2m+3 /t3/2g2m+3 

where 

2 

In this double scaling limit the solution u{x,t,e) of the KdV equation il.l(]\} admits the 
asymptotic expansion 

2 2 / _2m + 2 2m \ 4 

u{x,t,e) = Uc--e'2"^+sqU 2m+3 7-o(x, t,e),e 2m+3ri(i, e),0,0, . . . , Oj +C'(e2'"+3), (1.21) 

where 

X — Xc — Qucit — tc) , , 3(t — tc) , , 

ro{x,t,e) = ^—^ ^, n{t,e) = -^^. (1.22) 

Here q{s,ti,t2, ■ ■ ■ ,tm-i) is a solution to equation \1.S\) which has properties (i)-(ii)-(iii) 
given in Theorem M.li 



Remark 1.6 In the case of a generic critical point where ?7i = 2, Theorem [L5] describes 
the main theorem proved in [5j. The result is new for m> 2. 

Remark 1.7 It is likely that the above result holds also for other equations in the univer- 
sality class of Dubrovin [TB], which contains among others the Camassa-Holm equation, 
the de-focusing nonlinear Schrodinger equation, and the KdV hierarchy. 



1.1.1 Correlation kernels in critical unitary random matrix ensembles 

The first conjecture about the existence of real pole-free solutions to the Pf" equations for 
even m was posed in random matrix theory [21 13] in the case where ti = . . . = tm-i = 0. 
However, it leads no doubt that also the general case with non-zero ij's is of interest when 
studying critical unitary random matrix ensembles. Consider the space of Hermitian nxn 
matrices with a probability distribution of the form 

— exp(-n try (M))(iM, (1.23) 

where y is a scalar real analytic function with sufficient growth at ±00, for example a 
polynomial of even degree with positive leading coefficient. The limiting mean eigenvalue 
distribution for random matrices in such an ensemble has the form [lOj 



d^v{x) = ipvix)dx, tpvix) = Jl yibj — x){x — aj)h{x), for x € L)j^i[aj,bj], 

i=i 

where h is a real analytic function. The number of intervals in the spectrum, the endpoints 
ttj, bj, and h{x) depend on the confining potential V{x) and can be found in terms of the 
unique solution to an equilibrium problem. Generically h does not vanish on the intervals 
[uj, bj], and in particular not at the endpoints aj, bj, so that the limiting mean eigenvalue 
density vanishes like a square root at the endpoints of the spectrum [36] . 

The two-point eigenvalue correlation kernel in the model (|1.23|) is given by 

g-f y(x)g-f 1/(2/) ^^_^ 

Kn{x,y) = {pn{x)pn^l{y) - Pn{y)Pn-l{x)) , 

x y Kji 

where the pj's are polynomials orthonormal with respect to the weight e~ on the real 
line; the leading coefficient oi pj is kj. Scaling limits of the two-point kernel give rise to 
well-known limiting kernels such as the sine kernel in the bulk of the spectrum (where 
V'y is positive) and the Airy kernel near an endpoint aj (or bj) where h{aj) 7^ (or 
h{bj) 7^ 0) [3 dU [12]. Near points in the spectrum where h{x) = 0, more complicated 
transcendental kernels appear in double scaling limits. Near singular edge points where 
h{bj) = 0,h'{bj = 0),h"{bj) / 0, a kernel related to the Pf equation was obtained [7j. 
Near higher order singular points where h{x) ~ c{x — bj)'^~^2 for m even (the case where m 
is odd cannot occur in unitary random matrix ensembles of the form (jl.23p ) , it is natural 
to expect a kernel related to the Pp equation. If one takes double scaling limits where 
V depends on n in a suitable way, the limiting kernel should depend on s and also on 
the m — 1 time variables ti, . . . ,tm-i- This is a difference compared to the situation in 
Theorem 11.51 where only one time variable ti is non-zero. 

Outline 

In Section [21 we prove the results stated in Theorem 1 1 . 1 1 ab out the equations in the first 
Painleve hierarchy. We will construct the pole-free solutions q in terms of a Riemann- 
Hilbert (RH) problem which depends on s,ti, . . . ,tm-i, and m. We will prove the solv- 
ability of this RH problem for m even and s,ti,. . . , tm-i G K. This will imply the absence 
of real poles for q. Using a Lax pair argument, we will explain the relation between the RH 

6 




Figure 1: The jump contour T. 



problem, the equations in the Painleve I hierarchy, and the KdV hierarchy. An asymptotic 
analysis of the RH problem will lead to asymptotics for g as s — >• ±00. 

In Section [3l we focus on Theorem II. 5i We will state a well-known RH problem which 
characterizes solutions to the Cauchy problem for the KdV equation, and rely on the 
techniques developed in [13l [HI [5] to transform this RH problem to an equivalent one 
suitable for asymptotic analysis as e — )• 0. The main new point is the construction of a 
local parametrix using the model RH problem associated to the pole-free solution to the 
P™ equation studied in the first part. 

2 Pole-free solutions to the Painleve I hierarchy 

In this section, we will construct the solutions q occurring in Theorem 11.11 in terms of a 
matrix RH problem. At several points, we will refer to [6] where the proof of Theorem ll.il 
has been given in the case m = 2, and where a more detailed exposition can be found. 
We consider the following RH problem. 

RH problem for ^: 

(a) ^ : C \ r — )• C^^^ is analytic, where F = |J)=i ^j is a contour consisting of four 
straight rays 

Fi : arg(C - zq) =0, F2 : arg(C - zq) = 0, 

F3 : arg(C - ^0) = tt, F4 : arg(C - ^o) = -0, 

each of them oriented from the left to the right, see Figure [TJ Here zq can be any 
real number, and 2m+3 ^ < 9 < it. 

(b) ^ has continuous boundary values ^I'±(C) when approaching (^ G F \ {zq} from the 



left or right, and they are related by the jump conditions 
'1 1 



^+(C) = ^-(C)io 1-, 



1 
1 1 



1 
-1 



CGTi, 


(2.1) 


C G Ts u r4, 


(2.2) 


CGTa. 


(2.3) 



*+(C) = ^-(C) 
^+(C) = ^-(C) 

(c) As C — >• cx), ^ has an expansion of the form 

(2.4) 
where h,q do not depend on (, where 



and 



. m— 1 . 

4 2m+3 v-^ 4 2j + l 



Remark 2.1 The jump contour F = T(zq, 6) of the RH problem depends on the values of 
zq and 0, but the RH solution corresponding to (zQ) Q) can be transformed directly to the 
RH solution for any other value of (zq, 6'), as long as 2m+3 '^ < 6' < it. This can be done 
by analytic continuation of the RH solution ^ across its jump contour. In this section, we 
will fix the values of zq = and 6 = f^^ifvr, but for the asymptotic analysis in Section [23| 
we will need to choose zq,0 more carefully. Near zq, we need to impose that ^ remains 
bounded in order to have a unique solution. 

The RH problem depends on s, ti, . . . , im-i- For fixed real values of those parameters, 
we will show that there exists a (unique) RH solution ^(C) aud constants h and q such 
that (|2.4p holds. We have ^ = ^(C; s,ii, • • • ,im-i), h = h{s,ti, . . . ,tm-i), and q = 
q{s,ti, . . . ,tm-i)- We will show that q satisfies the conditions stated in Theorem ll.il 

Remark 2.2 If m = and s = 0, this is a well-known RH problem which can be solved 
explicitly in terms of the Airy function and its derivative. For m = 2, C, = 60^'^C, 
T = —4 • 60~^' ti, and X = 60~^'^s, we obtain the RH problem studied in p] related to 
the pole-free solution to the fourth order equation (jl.Sp which is equivalent to equation 
(fOI) for m = 2. 



2.1 Solvability of the RH problem 

In order to prove the solvability of the RH problem for ^, we perform the transformation 



*(C) 



I i)*(«' 



e(C)cT3 



(2.7) 



where we have suppressed the dependence on s,ti, . . . ,tm-i in our notations. This leads 
to a slightly modified RH problem for $. 

RH problem for $: 

(a) <& is analytic in C \ F. 

(b) We have the jump relations 



$+(C) = <i>-(C) 
$+(C) = <i>-(C) 
$+(C) = <l>-(C) 



1 e-2^(^) 

1 
1 

em) 1 

1 
-1 



for C e Ti, 



for C G Ts U r4, 



for C G Tg. 



(2.^ 



(2.9) 



(2.10) 



(c) As C — )• oo, $ has an expansion of the form 

$(C) = (/+^ + o(r'))r^'^^iv. 

It is straightforward to verify by (j2.4p and (|2.7p that g is given by 



q = A 



1,12 



2A 



1,11- 



(2.11) 



(2.12) 



Following the general procedure developed in [2ll [231 E] > the solvability of a large class 
of RH problems is equivalent to the fact that a homogeneous version of the RH problem 
has no solution except the trivial zero solution. This follows from the observation that a 
singular integral operator associated to the RH problem is a Predholm operator of index 
zero. This procedure is explained in detail in |llt [28] and in [6] for the above RH problem 
in the case where m = 2. All of the arguments work for general m and general real values 
of the ij's and we will not repeat them here. In order to have solvability of the RH problem 
for <I>, it is sufficient to prove the following vanishing lemma. 

Lemma 2.3 (Vanishing lemma) Let s,ti, . . . ,tm-i G K, let ^q satisfy conditions (a) 
and (b) of the RH problem for $, and let in addition 



MO = o{c'/^) 



as C ^- oo. 



(2.13) 



Then $n = 0. 



Proof. The proof of the vanishing lemma is ahnost the same as in [Bl [TT] and goes as 
follows. We first collapse the jump contour T in the RH problem for $ to the real line. 
We do this by defining 

AiC) = MC)(l ~^), forO<argC<e, 



MO = '^oiO { ^9(0 J (J 0^), for0<argC<7r, 



M0 = MC)[_)eiO i), 



for — TT < argC < —9, 



A{C) = MC), for-^<argC<0. 

Then one verifies using the jump relations for <I>o that A is analytic in C \ M, and that it 
has the jump properties 

A+(C) = ^-(C)(^^2e_(C) j' forCeR-, (2.14) 

A+(C) = ^-(C)(^ ^ oj, forCeM+. (2.15) 

As C — >• oo, the behavior of A follows from the asymptotics for $o ^-^d the fact that the 
exponentials in the definition of A are uniformly bounded. We have 

A{C) = 0{C^/^), asC^oo, (2.16) 

uniformly for (^ E C \ M. Then we define Q{C) = A{(^)A^{Q, where A^ is the Hermitian 
conjugate of the matrix A. Q is clearly analytic in the upper half plane, because A is 
analytic in the upper and lower half plane, and it is continuous up to M. As ^ — >■ oo, we 
have Q{C) = 0{(~'^^'^), which implies that 

[Q+iOdC = (2.17) 

JR 

by Cauchy's theorem, and using (I2.14p - (12.15j) we obtain 

+ y^^A_(o(^' ^ qJ^?(o^c = o. (2.18) 

For real values of s, ti, . . . , tm-i, we have 



0+iC,s,ti,...,tm-i) = 9-{£,;s,ti,...,tm-i), for ^ < 0, 

10 



and adding (j2.18p to its Hermitian conjugate, we obtain 

Since e~^^'^''''*i'---'*™'-i) > 0, this implies that the first column of A_ is identically zero 
(because it is continuous). The jump conditions (|2.14p - (j2.15p can be used to prove that 
the second column of A+ vanishes as well. 

Writing out the jump relations for the entries of A for which we have not yet proved 
that they vanish, the RH problem for A decouples into two scalar RH problems. The same 
argument as in |1H Step 3 of Section 5.3] shows that those scalar RH problems have only 
the zero solution, and that A = 0. Consequently we have that $o = 0, which proves the 
vanishing lemma. □ 

As a consequence of the vanishing lemma, the RH problem for $, and thus also for ^ 
(one can invert the transformation defined by (j2.7p ). has a solution for real values of the 
parameters s,ti, . . . , tm-i- 

It follows from the general theory of RH problems that the subset of C™ of values 
(s, ti, . . . , tm-i) where the RH problem is solvable, is an open set, and this implies in 
particular that, for any {s,ti, . . . ,tm-i) £ W^, there exists a neighborhood in C™ such 
that the RH problem is solvable also in this neighborhood. Moreover if this neighborhood 
is chosen sufficiently small, condition (j2.4p is valid uniformly, see [6], and ^, q, and h 
depend analytically on each of the variables s,ti, . . . ,tm-i- Values of (s, ti, . . . , tm-i) 
where the RH problem is not solvable, correspond to poles of q. 

2.2 Relation between the RH problem, the Painleve I hierarchy, and 
the KdV hierarchy 



We will show that the function q appearing in the asymptotic expansion (|2.4p solves 
equation P™, and that it also solves the equations in the KdV hierarchy. This follows 
from Lax pair arguments which are rather standard, see e.g. |21l [22l [29] in general and 
[351 03] for the Lax pair corresponding to the first Painleve hierarchy. We recall the 
arguments briefly for the reader's convenience. 

2.2.1 Lax pair in s, (, and t^ 

Since ^ is differentiable in s, we can define 

L(C) = *.(C)*(C)~S ^(C) = ^c(C)^(C)"'- (2.20) 

The jump matrices for ^ do not depend on (^ and s, so L and A are entire functions in 
the complex plane, and because of the asymptotics ()2.4p . they are polynomials in ^. For 
L, one deduces directly from (j2.4p that 



where it was used in the latter equation that hs = q. This follows by developing ^^^^^ 
as ^ — >■ CO and imposing that the 12-entry of the ^~^-term is zero. For the ^-derivative, 

11 



we use (|2.4|) to conclude that ^ is a polynomial of degree m + 1, 

m+l 

A{C) = Y, M'^ (2-22) 

j=0 

where the ^j's depend on s,ti, . . . ,im-i but not on (". Let us take a closer look at the 
12-entry ^(C) := ^12 (C): by ([231) and (H^Q]), it has the form 

m—l 

/3(C) = /3(™+i)(C)+5^tfc/3(^)(C), (2.23) 

k=l 

where f3^^' is a polynomial of degree k — 1 which can be written in the form 

fc— 1 
/3(')(C) = E^^C^ (2.24) 

i=o 

where Cj is independent of the value of k in the above formula, and £_i = 4, £0 = — 4g. 
Similar formulas can be obtained for the other entries but are not needed. 
Now, as in |13], write 



-4(0 =(; i) (2.25) 

(it is easily verified that Trvl = 0). Then ^^^ = ^^^ implies the compatibility condition 

L(;-As + [L,A]=0, (2.26) 

and if we write this down entry-wise, we find 

Q, + /3(C + 2g)-7 = 0, (2.27) 

/3, + 2a = 0, (2.28) 

7, - 1 - 2a(C + 2q) = 0. (2.29) 

We solve the second equation for a and then the first for 7, and we substitute their values 
in the third equation. This gives 

^/?,,,-2/3,(C + 2(7)-2(7,/3 + l = 0. (2.30) 

The degree m + l term in ()2.30p is trivial, and the degree m term reads 4gs + £o,s = 0, 
which we knew already since Co = — 4g. From the term of degree j with 1 < j < m — 1, 
we get 

'Z'~'m—j — l,sss ^'~'m—j,s ^Q.'~"m.—j — l,s ^Qs'~'m—j — l — U, yZ.oi) 

and the constant term in ()2.30p gives 

m—l 



2g£'^_i - qsCm-i + Ytj[ -^^j~2,sss - 2g>Cj_2 - qs^j-2 j = 0, (2.32) 



^' A '~"m—l,sss 

3 
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or 

, m-l , 

l + —Cm+Y^tj—Cj.i = Q. (2.33) 

as ^-^ as 

Integrating this equation gives the Pj" equation (ll.2p . One shows using the asymptotic 
behavior of solutions to the Schrodinger equation ^^ = L^ that the constants of inte- 
gration when integrating (j2.33p and (I2.3ip are zero. This proves that q solves the P™ 
equation (|1.2p . 

The jump matrices for ^ are also independent of ti, . . . , im-i, and consequently B^^> = 
^ij,^^^ is a polynomial of degree fc + 1 in C. Exploiting the compatibility of the i^- 
derivative with the s-derivative, an analogous argument as before leads to the time flow 

2.3 Asymptotics for q 

We will now analyze the RH problem for ^ asymptotically as s — )• itcxD. We will use the 
Deift/Zhou steepest descent method to obtain asymptotics for \I'. For m = 2, this analysis 
has been done in [6], see also |301l31j. For m, > 2, the general approach remains the same, 
but in particular the construction of the 5- function and the contour deformation are more 
delicate. 

2.3.1 Re-scaling of ^ 

Until now, we have always considered the RH problem for ^ corresponding to 20 = and 
6 = 2m+3 '^' ^^^ Remark 12. 1[ For the asymptotic analysis of the RH problem, we will need 
other values of zq and 6 which we will specify later. 
Define 

^^'^^={1 o)^(NI^C;5,il,•••,tn^-l). (2.35) 

RH problem for Y: 

(a) Y is analytic in C \ r(zo,^), where zq = \s\ "^+'^zq. 

(b) We have the jump conditions 

y+(C) = r_(C)(j I), forCGTi, (2.36) 

y+(C) = r_(c)(| °), forCGr2ur4, (2.37) 

y+(C) = F-(C)f° iV for C era. (2.38) 



T 0/ ' 
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(c) As C — )• oo, we have 



2m + 3 . 



(2.39) 
where 



2m+3 



0(C;S,tl,...,tm-l) = |s| 2m+2 0(|s|,n+l^;s,t^,...,t^_^) (2.40) 

. m— 1 , 

2m+3 , , I/O ■r — -\ 4 , , j-m-1 2j + l 

C^ +sgn(s)C'/2 + ^ irr—rtj\s\ -+1 C^- (2.41) 

i=i ^^ + ^ 



2m + 3 

The matrix Ai in (jTMj) is the same as in (plT]) . so (fZT2]) holds. 

2.3.2 Construction of the (^i-function 

We proceed with our analysis in the case where ti, . . . , tm_i = 0. Straightforward modifi- 
cations described in Remark 12.41 allow us to treat also the general case where the tj's are 
fixed. We search a g'-function g = g{(^; s) of the form 

5(C) = (C-^o)'/'p(^), (2.42) 

where p{z) is a polynomial of degree m. Its coefficients and zq are uniquely determined 
by the condition 

5(C) = ^(C) + 0(r'/'), asC^oo. (2.43) 

We have 

^(C) = ir^C"^ + sgn(.)C^/^ (2.44) 

2m + 3 

and because 

(C - zo)-'/' = C"' h + f; ^^It^-^C-M , as C ^ oo, (2.45) 

(I2.43P requires us to take p of the form 

The missing condition in order to have (|2.43p is 
, , /2'"-i(m + l)!\^ 
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With this choice of g, we have 

2m + 3 . °° 

e\s\^^(9(O-0(Oh^ = I + Y,dk<ylC-^'\ as C ^ oo, (2.48) 

fc=i 

where the coefficients dk can be calculated but are not important. Since the determinant 
of the left hand side of (j2.48p is 1, we have 

d2 = \dl (2.49) 

Remark 2.4 If the tj's do not vanish, we need to modify p and zq, but (j2.43p still deter- 
mines the coefficients Cj of the polynomial 

4 "^ 

P(-) = ^;^^^-o™E^^-^""' (2-50) 

uniquely. We don't need their explicit form, it is sufficient to know that 
(2j + l)!! 



2 



2Kj\ 



+ 0{\s\ ^+1), as|s|-^oo, (2.51) 



zo = -sgn(s) ^ /• +0(|sr^^), as|s|^oo. (2.52) 

\ [Im + 1)!! / 

To prove Theorem 11.11 we only deal with the case where ti = . . . = tm-i = 0, but the 
entire analysis done below can be easily generalized as long as the tj's remain bounded. 
Formula pT9|) wih then follow from (fTTI]) together with ([232]) . 

Proposition 2.5 We have 

giC) > 0, for C > zo, (2.53) 

Im5V(C)>0, forC<zo. (2.54) 

Proof. For the first equality, we ffist prove that p{z) > for z £ (— oo, — 3/2]U [— 1, +oo). 
Since the coefficients of p are positive, this is clear for z positive. For — 1 < z < 0, we have 

p{z) = — Z^ ({Cm - Cm-l\z\) + Z^ {Cm-2 - Cm-skl) + • • • + z"^'^ {C2 - Ci\z\) + CqZ™) , 

2m + 3 

and all the terms in this expression are positive since the Cj's increase with j. For z < 
—3/2, a similar argument shows that p{z) > (the terms in the sum are alternating 
and their absolute value increases with the degree). This implies that g{C) > for ^ G 
(— oo,— 3/2] U [— l,+oo]. In the case where zq > {or s < 0), this implies (j2.53p . for 
zq < 0, we still need an estimate for -f- G (-3/2,-1), or ^ G (— zq, — f-zo)- We already 
know that g is positive at both endpoints of the interval, so it suffices to prove that g(() 
is monotonic on {—zq,—^zq). 

Write <7'(C) = (C-^o)'/'g(|;) with 

(2j-l)!! 



q{z) = 2z^[z^ + ^bjz^-^], b 



J 



2Jjl 



15 



We will show that q{z) is strictly positive for all ^; G M, which implies (j2.54p . and also 
the fact that g{(^) is monotonic on {—zq,—^zq), which completes the proof of the first 
inequality. 

For z G (—00, — 1] U [— ^, +00), it is not difficult to see that q{z) is positive. Indeed, the 
terms in the sum are either positive or alternating and monotonic with increasing degree. 



For — 1 < 2; < — i, write 



z, 



^2-m 



z q[z) 



m 

+ r + l)+iZ (^^.-i^-^'^^^ + b,,z-^^+^) . (2.55) 



The first part is bigger than -j^, and each part of the sum reaches its minimal value on 



[—1,-1/2] at —1. Since the 6j's are decreasing, rearranging the alternating terms, we 
obtain the estimate 



z^-'^qiz) >T7: + i2 (-^2fc-i + &2fc) >—- — +bm>0. 



n) ^2-m 



16 



fc=2 



16 16 



n 



The above result enables us to apply the Cauchy-Riemann conditions, which leads to 
the following corollary. 

Corollary 2.6 There exists a 9 > such that Regf(C) < if arg(<^ — zq) = tt ±9. 

We still have the freedom to choose the values of 6 and zq that determine the jump contour 
T{zo,9) for Y. We take zq as in ()2.47p and 9 such that the above corollary holds. 



2.3.3 Normalization of the RH problem 

Define 

\(ii|s|2™+2 1 
so that we have 



n \ 2m+3 

^ 1 y(^)el^l^'^^9{CV3^ 



(2.56) 



RH problem for S: 

(a) S is analytic in C \ r(zo, 9). 

(b) ForCGr(zo,e), 

/ 2m+3 N 

1 e-2|^l^^^^s(C) 



S+(C) = S_(C) 



,0 



1 



/ 1 0^ 

5, (C) = S (C) \ 2m+3 



5+(C) = S_(C) 



1 
-1 



forCGFi, (2.57) 



for C G Ta U r4, (2.58) 



for C G Tg. (2.59) 
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(c) As C — )• oo, S behaves like 

5(C) = [/ + BiC^ + 0{C^)] \s\-^^^C^N, (2.60) 

and q is given by 

g = |s|^Si2^2-2|s|^5i,ii. (2.61) 

As s — )• ±00, the jump matrices for S tend to constant matrices except near zq: indeed on 
Fi, r2, and r4 we have exponentiaUy fast convergence to the identity matrix by Proposition 

2.3.4 Global parametrix 

In the limit s — )• ±00, if we ignore a small neighborhood of zq, the RH problem for S 
reduces to a RH problem with a jump only on (— oo,zo)- We can explicitly construct a 
solution P^°°> to this RH problem. 

RH problem for P^°°^: 

(a) p(°°) is analytic in C \ (—00, zq]. 

(b) We have 

pf)(C) = p(°°)(C)(_°^ J), forCG(-oo,Zo). (2.62) 

(c) As C — ^ 00, 

p(~)(C) = (I + 0{C^)) Isl'^^^^'C^N. (2.63) 

This RH problem can easily be solved explicitly: if we take 

P^^\(;) = \s\'^^{C-zoy^N, (2.64) 

one verifies that p(°°) satisfies the required conditions. The asymptotic condition (j2.63p 
can be specified as 



p(^)(C) = ^i+^a3 + 0{C')j \s\-^^^C^N. (2.65) 

We will show that P'°°) determines the leading order asymptotics of S and thus indirectly 
of the Pp solution q, by (j2.6ip . Therefore we first need to know that there exists a local 
parametrix near zq which matches with the global parametrix. 

2.3.5 Local parametrix near zq 

Let us fix a small neighborhood U oi zq, for example a small disk. Given m and sgn(s), 
we take U fixed for \s\ sufficiently large. Near zq, the ^-function vanishes like c(C — zqY''^. 
Following a well understood procedure, one can explicitly construct a local parametrix 
in U in terms of the Airy function and its derivative. We refer to for the explicit 
construction (and to HH [E] for similar constructions). The only thing that we need 
here, is the existence of a local parametrix which satisfies the RH problem 
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RH problem for P: 

(a) P is analytic in U \ T{zq, 9), 

(b) for ( E r(zo, 9)riU, P satisfies exactly the same jump conditions than S (see (I2.57P - 

dZSS])), 

(c) for C G dU, we have 

P{()p(^){()-^ = I + 0{\s\-^), ass ^±00. (2.66) 

2.3.6 Final transformation 

Define 

^^ \5(C)P(°°)(C)"\ forCGC\t/. ^ ^ 

Then R is analytic in the interior of U and across (—00, zq) because the jumps of S cancel 
against the jumps of the parametrices P and P'°°) . On the boundary of U, the matching 
of the local parametrix with the global parametrix, see (|2.66p . implies that R has a jump 
that is / + Odsl"^) as s — ;> ±00. 

RH problem for R: 

(a) R is analytic in C \ J^r, with S^j = (F \ [/) U dU, 

(b) R+iC) = R-{C)vr{0 for C e S^, where 

i;^((^) = / + 0(|sr^), as s ^ ±00, for C e 5f7, (2.68) 

t;^(^) = / + 0(e-"l"l(l^l+^)), as s ^ ±00, for C G S/j \ aC/, (2.69) 

(c) There exists a matrix i?i = Ri{s) such that i?(C) = I + RiC^^ + C(C^^) as C — )■ 00. 

It is a standard fact that the solution to a RH problem of this form (with small jump 
matrices and normalized at infinity) is close to the identity matrix |12j : we have 

R{C) = I + 0{\s\-'^), ass ^±00, (2.70) 

and for the residue matrix at infinity we have 

Ri{s) = 0{\s\-'^), ass ^±00. (2.71) 

Using (I2.65P and (I2.67p . one derives the identity 

Ri = Bi- ^(73, (2.72) 

and this implies that 

^1,11 = j + 0(|s|-i), 5i,i2 = 0(|s|-i), ass ^±00. (2.73) 

By (j2.6ip we obtain 

q{s) = — ^|s|™+i +0(151 ^+1), ass-;>±oo, (2.74) 

Now (II. 7p follows from (I2.52p . and the more general asymptotic formula (II. 9p follows from 
([232]) . 
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3 Critical behavior of solutions to the KdV equation 

In this section we will prove Theorem 11.51 and show that the pole-free solutions to the 
even members of the Painleve I hierarchy describe the critical behavior of solutions to the 
KdV equation. A RH procedure to obtain asymptotics for KdV solutions was developed 
in |13l [TJj . In [5] the method was used to prove Theorem 11.51 in the generic case where 
m = 2. For the sake of brevity and because many of the arguments are valid also for 
m > 2, we will refer to this paper at several points. In the RH analysis of the KdV RH 
problem, we will construct auxiliary matrix functions S, P^°°\ P, and R. They are not 
the same functions as in the previous section, we hope this does not cause any confusion. 

3.1 RH problem for the KdV equation 

Given initial data uo{x) satisfying the conditions specified in the introduction (i.e. uq(x) 
is negative, real analytic, has a single negative hump, and decays sufficiently fast at ±oo), 
we are interested in the solution u{x, t, e) to the Cauchy problem for the KdV equation 
(jl.lOp . The following RH problem characterizes u{x,t, e) at any time t > 0. 

RH problem for M: 

(a) M : C\M ^ C^^^ is analytic. 

(b) M has continuous boundary values M_|_(A) and M_(A) when approaching A G M\{0} 
from above and below, and 

/ 1 (\ \ 2ia(\;x,t)/e\ 

^+(^)=^^-(^)Ua;.K-(^-)A l-KA;e)P J' ^^ ^ < «' 

M+(A) = M_(A)c7i, cJi=(l J), forA>0, 

with a given by 

a(A; X, t) = 4t(-A)3/2 + x{-X)^/\ (3.1) 

The branches of (— A)^'^ and (— A)"*^'^ are analytic in C \ [0,+oo) and positive for 
A<0. 

(c) As A — )• oo, 

M(A) = (l + 0(A-i)) I |. (3.2) 



The solution M = M(A; x, t, e) depends on x, t, e. If r(A; e) is the reflection coefficient from 



the left for the Schrodinger equation e^-^/ + uo{x)f = \f with potential uq, then it is 



known that 



d 

u{x,t,e) = -2ie—- lim ( V-A[Mii(A; x, t, e) - 1] ) (3.3) 

ox A-i>oo ^ ' 
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Figure 2: The jump contour S5 after the transformation M >-^ S 



is the solution to the KdV equation with initial data uq{x) at time t > 0. Using certain 
smoothness and asymptotic (as e — )• 0) properties of the reflection coefficient, the RH 
problem for M can be transformed to a RH problem with modified jump matrices. We 
refer to [5] for the explicit construction of the function S which satisfies the RH problem 
stated below, with jumps on a deformed jump contour, see Figure [21 lenses are opened 
along an interval (—1 — S,Uc) for some small 5 > 0. The point Uc is the Hopf solution 
u{x, t) evaluated at the point x^ and time t^ of gradient catastrophe. 

RH problem for S: 

(a) S is analytic in C \ II5, 

(b) 5+ (A) = S-{X)vs{X) for A G S5, with 



vsW 



^1 m(A;6)eT<^W^ 
>0 1 



on Si, 

on S2 = Si, 

as A G (uc, 0), 



(3.4) 



(^m(A;6)e-T'^(^;-) ij' 

y iR{X;e) (1 - |r(A)|2)eT'^+(^;^)J ' 
(Ti, as A G (0, +00), 

and for A G (—00, —1 — 5), we have 

vs{X) = 1 + 0(e-t(l^l+^)), as e ^ 0, 0, 

uniformly in A for x,t sufficiently close to Xc,tc- 

(c) S{X) = {I + O(A-i)) ( .^ _.^ as A ^ 00. 

S can be expressed explicitly in terms of M [5], and it follows from this explicit expression 
that 



(3.5) 



u{x, t, e) = Uc- 2ie—Si^ii{x, t, e). 



(3.6) 
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where 

Su{X;x,t,e) = l + ^^^^^pd^ + 0{X-^), asA^oo. (3.7) 

The function k can be expressed in terms of the (analytic continuation of the) reflection 
coefficient and satisfies the important asymptotic property 

K{X;e) = l + 0{e), for A G Si U [n,0], as e ^ 0. (3.8) 

Furthermore (p depends explicitly on the initial data: 

^{X; X, t) = ^Jua -X{x -Xc- 6uc{t - Q) + 4(uc - Xf/^{t - Q 

+ pifiiC) + Gtc)Vl^Xd^. (3.9) 

If 



/f (^,) = /f (^,) = /f (n,) = . . . = ft\u,) = 0, /f +'Vc) / 0, (3.10) 

repeated integration by parts gives (since f']^{uc) + 6tc = 0) 

(/>(A;x,t) = Vuc - X{x -Xc- 6uc{t - Q) + 4(uc - Xf/'^{t - Q 

+ 7^^ r fi^^'^^OiC - xf-^dC. (3.11) 
(2m + 1)!! 7;, 

For any fixed neighborhood U of Uc, it was also proved in [5] that there exists 6 > such 
that 

jl + 0{e-'l), forAeS5\(^U(t.„+oo)), 
^5 A) = < , r ^ / nN\7/ as e ^ 0, (3.12) 

[tai + 0{e), foi Xe {uc,0)\U, 

if |x — Xc\ < 6 and |t — tc\ < 5. 

3.2 Construction of the global parametrix 

If we ignore the jump matrices that are small as e — t- and the jumps in a fixed sufficiently 
small neighborhood U of Uc, we obtain the following RH problem: 

RH problem for p(°°): 

(a) P(°°) : C \ [uc, +oo) -^ C2x2 is analytic, 

(b) P^°°> satisfies the jump conditions 

p(oo) ^ pM^^^ o^ (0,+oo), (3.13) 

pjoo) ^ .p|oo)^^^ on('Uc,0), (3.14) 

(c) p(°°) has the following behavior as A — )• cxd, 

1 1 

,i(-A)i/2 -i(-A)i 

This RH problem is solved by 

'1 1 



P(-)(A) = (/ + 0(A-1))(^, ^ /2 _,v_\u/2)- (3-15) 



p(°°)(A) = (-A)i/^(ue - A)-"='/^ r 



i — i 
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(3.16) 



3.3 Construction of the local parametrix 



We need to construct a local parametrix in a neighborhood lA of Uc- As e — )• 0, we have 
K,{\) = 1 + 0(e), and we will construct a function P, defined in lA^ which satisfies the same 
jump relations as S, but in the limiting case where k, is set to 1. 



RH problem for P: 

(a) P : Z7 \ S5 ^ C2><2 is analytic, 

(b) P satisfies the following jump condition onU r\ S5, 

P+(A) = P_(A)«p(A), 
with vp given by 



vp{X) 



.0 1 J' 

^•g-f<^(A;x,t) ^j 

i 0, 



as A £ El, 
as A G E2, 
as A G (uc;0), 



(3.17) 



(3.18) 



(c) in the double scaling limit where e — )• and simultaneously x — )• Xc, t — )• tc in such 
a way that 



X — Xc — QuJt — tc) —3(t — tc 
hm ,;::xo = To, lim — ,^ 



221+2, 1/, 



g2m+3/j3/2 

with A; given by (jl.20p . we have the matching 

P(A)P(°°)(A)"^ ^/, forAeaZY. 



Ti, ro,Ti e 



(3.19) 



(3.20) 



We will use the RH solution ^ = vj/C'") studied in Section [2] to construct the local 
parametrix P. First we transform the RH problem for ^ to a RH problem for <1> which 
models the jumps needed for P in an appropriate way. 

3.3.1 Modified model RH problem 

Define 



$(C;s,ii) = e"f^'^3^(C;s,ti,0,...,0)e 

for Im C > 0, and 

$(C;s,ti) = e-x-3^(^;s,i,,0,...,0)e' 



e(C;s,ti,0,...,0)<73 



-1 

1 



e 4 



0-3 



e(C;s,ti,0,-,0)a3gf ^3 



(3.21) 



(3.22) 
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for Im(^ < 0. We also write 

ms,h) = -^^i-cf-^ - ^hi-cf' + .(-c)^/^ (3.23) 

which is related to 9 in the case where t2 = ■ ■ ■ = im-i = 0, but with its branch cut on 
(0, +oo). One has the identities 

e = ie+, on(0,+cx)), 9 = ie, onTs, 6 = -iO, onT4. (3.24) 

Then it is straightforward to verify that <I> solves the RH problem 

RH problem for $: 

(a) $ is analytic for C G C \ f, with f = Ti U r2 n r4. 

(b) <^ satisfies the following jump relations on F, 

<I>+(C) = <^-(C)( ^ A, forCeri, (3.25) 

<I>+(C) = <&-(C)(J '" ^ 1, forcer^. (3.26) 

<^+{0 = ^AOi.J<,.,.,^, °1, forCGr4. (3.27) 



ze 



-2ie{(;s,ti) -^ j ' 



(c) $ has the following behavior at infinity, 



<I>(C) = ^(-C)-^'^^ ( \ J) (/ + ihas{-Cr'/' + OiC')) , (3.28) 



V2^ 



with the branches positive for C < and analytic off [0, +oo). 
We search for a parametrix P of the form 

2 2m+2 2m. 

P(A) = E(A;e)$(e~5^^^/(A);e~5^^^To(A;x,t),e ^^^^ri(A; t)), (3.29) 

where E, f, to,ti are analytic in ^. So we evaluate <&(C; s,ii) at the values 

2 2m+2 

^ = e 2;;r+3/(A), s = e~5;^^To(A;x, t), (3.30) 

2m 

ti = e 2;^^ri(A;t), (3.31) 

and we will construct /, tq, and ti in such a way that 

~ 2 2m+2 2m 1 

6i(e 2m+3/(A);e 2'"+3To(A;x,t),e 2m+3ri(A; t)) = -0(A;x, t). (3.32) 

This condition is satisfied if we define / by 

4 , , , , 2m+3 2™ /'"'= (m-l-^^ , s. , 2m+l 

/(A))^= ,, , ,,„ / ft^'\0{^-^)^d^, (3.33) 



2m + 3' ^ " (2m + 1)!! jx 
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n by 

-^n(A;t)(-/(A))i =4(t-y(^,-A)3/2, (3.34) 

and To by 

ro(A;x,t)(-/(A))5 = ^u^ - X{x - x^ - 6uc{t - t^)). (3.35) 

Indeed, summing (j3.33p - (|3.35p gives (|3.32p by p. lip . This defines /, ro,ri analytically 
near Uc, and we have 

2 

/(ue) = 0, f'iu,) = (^-^__^/f+i)(n,)J = fe > 0, (3.36) 

riin.) = -%/^, (3.37) 

X — Xc — 6Uc(i — ic) 

ro(nc) = ^jjy^ . (3.38) 

Since / is a conformal mapping from a neighborhood of Uc to a neighborhood of 0, we can 
choose the lenses of the jump contour for S in such a way that /(S5 OU) C F. Then for 
any analytic function E near Uc, P satisfies the required jump conditions on Sg n Z// (see 
(|3.18p ). but we also need the matching (|3.20p . which has to be valid in the double scaling 
limit where e — )• 0, x — t- Xc, t — t- tc in such a way that (I3.19P holds, or in other words 

2m+2 2m 

lime 2m+3 7-Q(tt^; x,i) = To, lime '^^+'iTi{uc;t) = ti. 

2m+2 2m 

If ZY is sufficiently small, (e '^"^+^TQ{\;x^t)^e '^'^+'-^Ti{\]t)) will lie in a small complex 
neighborhood of (ro,Ti) for A G dlA. By ()3.20p . we have 

P(A)P(-)(A)-1 = -Li^(A)(-e-2^/(A))-i'^3 (^_1^ 1^ 

('/ + i/i(j3e^^(-/(A))-^/2 + 0(era)')p(°°)(A)-\ (3.39) 

2m+2 2m 

as e —)• 0, where /i = /i(e ^■^+^TQ{\;x,t),e '^'^+^Ti{\]t)). If we define 



i5;(A) = i=p(-)(A)(} ^^)(-e- 



2m+. 



3/(A))2-3, (3.40) 



it is easily verified that E is analytic in lA and that we have 

P(A)P(°°)(A)~^ = P(°°)(A) 

(/ + i/ia3e^^(-/(A))-^/2 + c'(gOT)') p(°°)(A)-\ (3.41) 

in the double scaling limit, for A G dU. 
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for A G S5 \ Z^, 


(3.44) 


for XgJ^sDU, 


(3.45) 


for A e dU, 


(3.46) 



3.4 Final RH problem 

1 
We define R in such a way that it has jumps that are uniformly I -\-0{e 2'"+3 ) in the double 

scaling limit: we let 

[S{X]x,t,e)P{X;x,t,e)~\ asXeU. 
Then, using the fact that 
vs{X;x,t,e)vp (A;x, t, e) = / + 0(e), uniformly for A E Z^ n S5 as e —)• 0, (3.43) 

one can verify that R solves a RH problem of the following form. 

RH problem for R: 

(a) R is analytic in C \ (E5 U dU). 

(b) R has the jump condition i?_|_(A;x, t,e) = R^{X;x,t,e)v£i{X; x,t, e) for A G II5 U OZY, 
where 

VR{X;x,t,e) =/ + 0(e-f), 
'(;K(A;x,t,e) = / + 0(e), 

VR{X;x,t,e) =/ + 0(e5^), 
in the double scaling limit where e — )• 0, x — )• a^c, t ^ tc and simultaneously 

2m+2 2m 

e 2m+3 7-o(nc;x, t) — )• To and e 2m+3-7-i(ttc; t) — )■ ri. 

(c) As A — )• 00, we have 

i?(A;x,t,e)=I + :^^i%^ + 0(A-2). (3.47) 

A 

On dli with clockwise orientation, the jump matrix has the form 

VR{X]x,t,e) = / + t'i(A;a;,t)e2^ + 0(52^^^), (3.48) 

with 

v^{X-x,t) = ih.{-f{X))-^'^P^^\X)a;P^^\X)-\ (3.49) 

This is a meromorphic function in li with a simple pole at Uc, the residue is given by 

_2m+2 2m _, ,„ /Q 2^\ 

Res(vi;nc) = -/i(e 2™+3To(uc;a;, t), e 2m+3ri(uc; t))A: ^/^ ( ^j . (3.50) 

Then as in [SJ Section 4] one can conclude that 

Ri{x,t,e) = e2^^^Res(t;i;'Uc) + C'(e2^^^), (3.51) 
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and by (j3.6|) and (J3.38P this leads to 

d 

u{x,t,e) = Uc-2e—Ri^i2{x,t,e) 

1 /o 2 dTniUr' X t) 2m + 2 '2m 4 

= Uc + 2k~^'^e^^^^ \ ^ g(e 2;^^ro(nc;x, t), e 2^^^ri(uc; t),0, . . . ,0) + Ofez^?^^ 

ax 

_, 2 2m+2 2 m 4 

= Uc-2k e2™+3g(e 2m+3ro(nc;x, t), e 2m+3ri(uc; t),0, . . . ,0) + ©(ez^+a), 
which proves Theorem 11.51 
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